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Model ol angular bremsstrahlung of photons emitted during a-decay is presented. A special em- 
phasis is given on development of unified formalism of matrix elements in the dipole and multipolar 
approaches. A probability of the emission of photons calculated on the basis of the multipole model 
without any normalization on experimental data (i. e. in absolute scale) is found at 90° of the angle 
i9ci7 between directions of motion of the a-particle (with its tunneling under barrier) and emission 
of photons to be in a good agreement with the newest experimental data for the Po, Po, and 
^^'^Ra nuclei. The spectrum for ^"^^Cm is found at i9ct7 = 25° to be in satisfactory agreement with 
high limit of errors of experimental data of Japanese group. A comparative analysis for the spectra 
calculated for ^^"Po by the multipole and dipole approaches in the absolute scale and with normal- 
ization on experimental data is performed. The emission of photons from the internal well in the 
dipole approach is found to be not small while the multipolar approach does not show such a strong 
dependence. Distribution of the bremsstrahlung probability on the numbers of protons and nucleons 
of the a-decaying nucleus in selected region close to '^^"Po is obtained. An unified formula of the 
bremsstrahlung probability during the a-decay of arbitrary nucleus expressed directly through the 
Qct-value and numbers Ap, Zp of nucleons and protons of this nucleus is proposed. 
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I. INTRODUCTION 



Last two decades many experimental and theoretical efforts have been made in investigation of nature of 
bremsstrahlung emission which accompanies a-decay of heavy nuclei (see references in dll). A key idea of such 
research consists in search of new information about dynamics of the a-decay (and dynamics of tunnehng) which is 
supposed to be extracted from analysis of the measured bremsstrahlung spectra. Tunneling times in nuclear collisions 
and decays are extremely small, close to nuclear time. This fact results in impossibility to test non-stationary methods 
of tunneling experimentally. However, researchers open new ways how to obtain new information about the dynamics 
of nuclear processes. An increasing interest in study of the bremsstrahlung processes accompanying the a-decay could 
be explained by this idea mainly: to find a new approach how through analysis of the bremsstrahlung spectra to 
"measure" dynamics of the a-decay (perhaps, in its first stage), to estimate duration of tunneling of the a-particle 
through barrier. 

Many approaches for description of the bremsstrahlung emission during the a-decay have already been developed 
where models with semiclassical spherically symmetric description of the a-decay are prevailing (see [1, [3, H, , also 
calculations of the spectra see in [7|). In comparison with fully quantum approach, the semiclassical one allows to work 
with characteristics and parameters, physical sense of which is natural that simplifies this task, allowing to understand 
studied questions easier. Enough well description of experimental data has already been achieved in such approach, 
where one can note a resent success in agreement between theory and experiment for the controversial nucleus ^^°Po 
M, . Perspectives are certain in study of dynamics of the a-decay with some analysis of the bremsstrahlung 0, 
[ll|, in study of dynamics of tunneling in the a-decay ^H, [3, [fl ITsl. in research of pecufiarities of the polarized 
bremsstrahlung during a-decay and influence of electron shells on it [l6j|. in effect called as Munchhausen effect 
which increases penetrability of the barrier due to charged-particle emission during its tunneling and which could 
be interesting for further study of the photon bremsstrahlung during tunneling in the a-decay. However, the fully 
quantum approach (starting from [Tsj and then [lol . [20j ) seems to be the most accurate and motivated from the 
physical point of view in description of emission of photons, to be the richest in study of quantum properties and new 
effects. Among the fully quantum approaches a model proposed for the first time by Papenbrock and Bertsch in [l9j 
has been developing the most intensively, where wave function of photons is used in the dipole approximation. In such 
dipole approach the matrix element is calculated with higher convergence and without visible decrease of accuracy, 
that makes this problem to be studied for many researchers in the fully quantum consideration. 

If intensity of bremsstrahlung was varied enough visibly at change of the angle, then the emission of photons would 
take influence on d yna mics of the a-decay essentially and, therefore, change some its characteristic. From such point 
of view discussions [2ll |22| | show a way for obtaining the new information about the a-decay: through angular analysis 
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of the hremsstrahlung during the a-decay. But for such researches a model of description of the bremsstrahlung in 
the a-decay which takes a value of the angle between the directions of the a-particle propagation (or tunneling) and 
the photon emission into account, should be constructed. In such direction two independent groups have developed 
the angular formalism in the semiclassical dipole approach including quadrupole term and in the fully quantum 
approach [Ij, i23l . [2^ . [25j where wave function of photons was mainly expanded by spherical waves. In particular, 
the good agreement for the ^^"Po nucleus between experimental data and the calculated spectra in the approach [1, 0] 
was obtained, neglecting by emission of photons from the internal nuclear region before the barrier. From such results 
it could follow that such internal emission of photons is extremely small and does not influence on the spectra. In 
such a case, it looks to be impossible to extract any useful information about processes of a-decay inside the nuclear 
region before the barrier from experimental data of the bremsstrahlung. But, in order to clarify this the emission of 
photons from internal region should be estimated on the basis of fully quantum calculations based on the realistic 
potential between a-particle and daughter nucleus. A question of influence of emission from the internal region of the 
total bremsstrahlung spectra on the basis of realistic a- nucleus potential has not been studied yet. 

Taking into account expansion in spherical waves for description of angular correction of the wave function of 
photons f23j and realistic form of interaction between the a-particle and the daughter nucleus [H, we achieved a 
little better agreement between such fully quantum approach and later obtained experimental data [8| for the ^^"^Po 
nucleus (for explanation see Fig. 3 in [24| and discussions here) . Results in descriptions of the newest experimental 
data [1, 0, UHl for the ^^''Po and ^^^Ra nuclei seem to be enough good also, where we have been achieving agreement 
between theory and experiment up to 765 keV. But, a multipole approach (started from [20j . then in H^) seems to be 
the most accurate and corrected in spatial description of the emission of photons during the a-decay. Such approach 
turns out to calculate probability without any normalization relatively experimental data (i. e. in co-called absolute 
scale), and achieves enough good agreement with them. This peculiarity adds power of prediction to the multipolar 
approach: it allows to study the bremsstrahlung during the a-decays of arbitrary nucleus, to estimate the emission 
from the internal region, to predict new spectra. 

In ^23| dependence of the bremsstrahlung probability on the electrical charge of the daughter nucleus was analyzed. 
But, it is unclear how much energy of the a-particle takes influence on the photons emission. Calculating the 
bremsstrahlung probabilities for the different isotopes with different mass numbers, I have seen that this probability 
is determined by combination between numbers of protons and neutrons in the nucleus under decay, rather than 
by the electric charge of this nucleus, and one come to conclusion about direct dependence of the bremsstrahlung 
probability on the effective charge Z^s- It could be interesting to clarify whether any other parameter or characteristic 
exists, influence of which on the probability is essential. More intriguing task has been appeared: to construct an 
unified formula of the bremsstrahlung probability during the a-decay of the arbitrary nucleus which is directly expressed 
through all these parameters and characteristics. But, whether is it possible to describe the bremsstrahlung spectra for 
all different nuclei by only one formula in general? Whether is it possible to describe the bremsstrahlung spectrum for 
only one arbitrary nucleus with very high accuracy inside the energy region of the photon emitted used in experiments ? 
To answer on such questions, it should be desirable to use the model which is the most accurate and corrected in 
description of this reaction. 

This paper answers on these questions and it is organized so. In Sec. II the model of the bremsstrahlung accom- 
panying the a-decay is presented, where emphasis is made on the angular formalism in framework of the dipole and 
multipolar approaches of the matrix elements and calculation of the bremsstrahlung probability. In Sec. IIII Al the 
model is tested on experimental data 24, 25] for the ^^''Po and ^^^Ra nuclei. The spectrum for ^^''Cm is added 
to such a picture, compared with high limit of errors of experimental data [tI, [1^. In Sec. IIII Bl distribution of the 
bremsstrahlung probability on the numbers of protons and nucleons of the a-decaying nucleus in selected region close 
to ^^°Po is presented. In Sec. IIII Cl the angular spectra for the ^^"Po nucleus by the multipolar model are presented. In 
Sec. IIII Dj the formula of the bremsstrahlung probability during the a-decay, based only on the Qa-value and numbers 
of protons and neutrons of the decaying nucleus, has been constructed. Inside region of the a-active nuclei from ^"^^Te 
up to the nucleus with numbers of nucleons and protons Ap — 266 and Zp = 109 with energy of the photons emitted 
from 50 keV up to 900 keV good agreement has been achieved between the spectra, obtained on the basis of the 
multipole model, and the bremsstrahlung spectra obtained on the basis of the proposed formula. At finishing, results 
are summarized. 



II. MODEL 

A. Matrix element of emission 

We define probability of the bremsstrahlung emission during a-decay of nucleus in terms of transition matrix 
element of the composite system (a-particle and daughter nucleus) from its state before emission of photon (called 
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as initial i-state) into its state after such emission (called as final f -state). In this paper, I shall use the definition of 
the matrix element like (2.11) in [26 | (in the first correction of the non-stationary perturbation theory with stationary 
limits to — —GO and ti = +00, and with normalization \C\ 1): 



(1) 



where 



P{kt,kf) 
pih, kf) 



(2) 



-zkr 



dr 



k^^ j rf{r) e-'^^'-^Mr) dr 



and V'i(r) = l^i) and V'/(r) = are stationary wave functions of the decaying system in the initial i-state and final 
/-state which do not contain number of photons emitted, Zcs and m are effective charge and reduced mass of this 
system, e*^") are unit vectors of polarization of the photon emitted, k is wave vector of the photon and w = k = |k|. 
Vectors e*^"-' are perpendicular to k in Coulomb calibration. We have two independent polarizations e^^^ and e*^^^ 
for the photon with impulse k (a = 1,2). One can develop formalism simpler in the system of units where fi — 1 
and c = 1, but we shall write constants U and c explicitly. Let's find also square of the matrix element aji used in 
definition of probability of transition. Using the formula of power reduction of 6 -function (see [29j . § 21, p. 169): 



[d{w)]^ = 5{w) S{0) = S{w) (27r)-i / dt = d{w) {27r)-^ T, 



we find (T — > -|-oo is higher time limit) : 

that looks like (4.21) in ^ (with accuracy up to factor (27r)2) and like (42.5) in [s^ (exactly, see § 42, p. 189). 



(3) 



(4) 



B. Linear and circular polarizations of the photon emitted 



Rewrite vectors of linear polarization e^"^ through vectors of circular polarization with opposite directions of 
rotation (see ref. [HI, (2.39), p. 42): 



Then p{ki,kf) can be rewritten so: 
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(5) 



(6) 



l±i 



\ 2, a=l,2 



(7) 



C. Different expansions of the vector potential A 

In order to simplify calculations of the matrix element p{ki,kf) with taking into account angular correlation between 
vectors k and r, function e"*"^"" connected with the vector potential A of the electro-magnetic field of the daughter 
nucleus should be expanded into basis of functions. At present, only three different types of expansion have been used 
in this problem: 

• The expansion in the dipole approximation: 

e*"" = 1 + « kr + . . . (8) 

where the first item is found to be mainly used correctly in further calculations of the matrix element; 
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The expansion by the spherical waves (in direction of papers [U, 0, [H, [U [2^ ) : 

g^kr^'^ ./(2; + i)P;(cos0„^)j,(fcr), (9) 



where Oa-y is angle between vectors k and r (kr — kr cos 
• The expansion by electric and magnetic multiples (see ref.^J, (2.106) in p. 58): 



£,^e'''^ = ^l^/2^Y.^2l + Dl^{^,9,Q) ■ [A,,(r, M) + A,,(r, S)j , (10) 



where (see ref.flU, (2-73) in p. 49, (2.80) in p. 51) 
Ai^{v,M) = ji{kr) Tu^^{nph), 



^iu{v,E) = yiJ-lj,_i(fcr)T,,_i,^(np;,)- y^^^ ^^^^ 

Here, Kiu{v,M) and Ai,j(r,E) are magnetic and electric multipoles, ji{kr) are spherical Bess el functions of order I, 
Tiii^^{n) are vector spherical harmonics, 9i, 62, O3 are angles defining direction of vector k relatively axis z in selected 
frame system. 

In this paper, I shall use the multipole expansion assuming that it gives us the most reliable approach to describe 
the angular emission of photon during the a-decay. Matrix- function Dl;*{(p, 9, 0) defines direction of vector k relatively 
axis z in the frame system for r: angles ip and 9 point to direction of vector k, but not the vector r. The functions 
T;i',i/(n) have the following form (^0 — Oj see ref.l3l|, p. 45): 

Tj;,m(n)= ^ (Z, 1, j I TO - ^,/i,m) r/,„_^(n) (12) 

Ai=±l 

where (?, 1, j | to — /i, /i, to) are Clebsh-Gordon coefficients and Yim{9,(p) are spherical functions defined, according 
to [13] (see p. 119, (28,7)-(28,8)). 



D. Approximation of the spherically symmetric Q-decay 

Now we shall study the a-decay in the spherically symmetric approximation. We orientate the frame system so 
that axis z will be parallel to the vector k and 

Dl^{^,9,0)=S^,. (13) 

In the spherically symmetric approximation, wave functions of the decaying system in the initial and final states are 
separated into the radial and angular components, and these states are characterized by quantum numbers I and to. 
We shall be interesting in such photon emission when the system transits to superposition of all possible final states 
with different values of the magnetic numbers toj at the same orbital number If. Let's assume that in the initial 
state we have h = mi = and the radial component of wave function (pf{r) does not depend on mf for selected //. 
We write wave functions so: 

tpi{r) = ipi{r)Yoo{nl.), 

^f{v) = vf{r) Y.Yi^.m{nl) (14) 

m 

and obtain: 

p {h, kf)^V2^Y^ {-lY V21 + 1 
I 

where 



Pf 



Wi 



(15) 
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and 

+00 



< = I drj dnrH}{v) (|:V^.(r)) A,;(r,M), 



+ 00 

Pf, = j drj dQr^rfir) (|:^.«) A,;(r,i?). 


Using gradient formula (see (2.56), p. 46 in (3l|): 



a , , , , I fdf l + l \ , ^ l + l fdf I , 



we obtain: 



, . ^ d(p,{r) 



and then we calculate the matrix components: 

Pl \ = — lM{lf,lphJph) ■ J{l,l), 



-r lEilfJphJph — ^) ■ J{lf,lph — ^) + \ r,, 1 lEilf,lphJph + i)-J{lfJph + i), ^ ^ 



where 



Jilf,n) = I ^*{l,r)^^Ukr)r^dr, 



lM{lf,lph,n) = I i^;,„(n;)Toi,o(n;)T;;^„^^(np;,)dr!, 

lE{lfJph,n) - f'f^ [^hrn'^''r)ToiAK)T^l,nA^ph)dn. 

LI J 



Ai=±l 

Using the following value of the Clebsh-Gordon coefficient (see Appendix A) 



from (in]) and IH]) we obtain: 



/i— — 1,1 



(17) 



V'.(r) = -^^Toi,o(n;), (19) 



(21) 



(110|l,-l,0) = (110|-l,l,0) = yi (22) 



(23) 



and for the angular integrals for transition into the superposition of all possible final /-states with different at the 
same If from eq. (|2ip we obtain: 



hi{lfjph,n) = E E ("'l'U|M^Ai',Ai',Ai) y yi*m(n^)^i,-M'«)'5C,M-M'("p'') 

/i — ±1 /i' — ±1 

lE{lf,lph,n) = E E (">1, V|M-Ai',Ai',Ai) y" ^;*m(n^)^i,-M'«)^n,M-M'("p'') '^^^ 
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E. Calculations of the angular integrals 



Let us analyze a physical sense of vectors nj,, n(. and Uph- According to definition of wave functions and 
'0/(r), the vectors n* and determine orientation of radius- vector r from the center of frame system to point where 
this wave functions describes the particle before and after the emission of photon. Such description of the particle 
has a probabilistic sense and is fulfilled over whole space. Change of direction of motion (or tunneling) of the particle 
in result of the photon emission can be characterized by change of quantum numbers I and m in the angular wave 
function: Foo(nr) Yimin() (which changes the probability of appearance of this particle along different directions, 
and angular asymmetry is appeared) . The vector Uph determines orientation of radius- vector r from the center of the 
frame system to point where wave function of photon describes its "appearance" . Using such a logic, we have: 

Uph ^ nl = n{. = Ur. (25) 

As we use the frame system where axis z is parallel to vector k of the photon emission, then dependent on r integrant 
function in the matrix element represents amplitude (its square is probability) of appearance of the particle at point 
r after emission of photon, if this photon has emitted along axis z. Then angle 9 (of vector rir) is the angle between 
direction of the particle motion (with possible tunneling) and direction of the photon emission. 
Let us consider the angular integral in ((24)) over dfi. Using ([25)) . we find: 



, .w+„_^'+i+^ / 3(2^ + l)(2n+l) (^-1)! (n-|m + //|)! 

^"^^ ' V 32^ (ITT)! (n+|m + M'|)! ^ ^ 

Pi (cos 6) Pi (cos 6) p1™+^' I (cos 9) ■ sin 6d9dLp, 

where P[^{cos9) are associated Legandre's polynomial (see [30|, p. 752-754, (c,l)-(c,4); also see [3l| (2.6), p. 34) and 
the following restrictions on possible values of m and If have been obtained: 

m = = ±1, lf>l, ^^ > - m'I = l"i + A^'l- (27) 

Let's introduce the following differential matrix elements dpf^ and dpf dependent on the angle 9 : 

dpf 



sin 9 d9 

dpf 
sin 9 d9 



m=±l fi' = ±l 




'--^J{lfjph-l) J2 ^hJpH.hu-~l fljJpH-l^^) 



(28) 



m— ±1 /J.' — ±1 



where 



CT.Cn - (-l)'^+"+^-^'+^(n,l,;,,|-m-M',M',-Hx 



X 



' (2^/ + l)(2n-|-l) {lf~l)\ {n-\m + ti'\)\ 
327r (// + !)! (n+ |m + /i'|)!' 

ir.niO) - Pl^{cos9) Pl{cos9) Pt^^'\cos9). 



(29) 



One can see that integration of functions (pS)) by angle 9 with limits from to tt gives the total matrix elements pf 
and pf exactly for transition into superposition of all possible final states with different my at the same I / . 

We shall find the matrix element at the first values of If and Iph- We have If = 1, Iph — \- Calculating coefficients 

C{"^ and functions (9) (see Appendix IB] and IC]) . from eq. (|^ we obtain: 



"^^i' - ^ -1^ J{l,l)-sin^9cos9, 



sin9d9 8 V tt 

dpf _ ,1 ./2 , .1/1 



sin 9d9 8 



• J(l, 0) • sin^ 9 + i^^- J(l, 2) • sin^ ^ (l - 3sin2 



(30) 



Integrating these expressions over angle 9, we find the integral matrix elements: 

= 0, = i 



(31) 



F. Dipole approximation and Fermi's Golden rule 

1. Vector potential A and integrated matrix element 
The matrix element p (ki, kf) in the dipole approximation of the vector potential A has a form: 
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(32) 



Applying the following transformation at li ~ introduced in [l^ (see new analog in a general case in Appendix IE)) : 



d 






dVir) 




dv 


1 Wf, 


dr 





we obtain: 



p{k„kf) = — J2 KC^ I rf{r)Mr) 



dV{r) 
dr 



dr. 



(33) 



(34) 



In the spherically symmetric approximation of the a-decay we use the wave functions in form (|14p and then the matrix 
element transforms into 



ih, kf) = — J2 E '^mC / ^''^^ I ■ ^/W yi]mMi) ■ vr{r) ■ 



dVjr) 
dr 



(35) 



After use of the gradient formula (fT9|) at = and taking into account form ([23|) for the spherical function Toi,o(nr), 
this matrix element is separated into radial and angular integrals: 



p{ki, k 



f) 



1 11 

Wf 



or (f±i and ^^i ^-re orthogonal vectors) 
where 



-\/| E '^A. E C^M' / ^}{r)^.{r)^^r'dr f Y*„^^inf) Y,,^,,in^) dQ (36) 

* u=-l.l u' = ±l i ^ rnt 



p {ki,kf) = Jdip (//) • /dip {If), 



(37) 



+ 00 



Jdip (If) / ^*f{lf,r) ip-,{r) 



dVjr) 2 , 
— r dr, 

dr 



^dip(//) = ^ 5] ^/i^y'y,;„,^(n;)yi,_4n,)da 



(38) 



After use of formula ([25]) for vectors and the angular integral /dip {If) is non-zero only at 

TO/ — — ^ = ±1, If — 1 

and equals to 

/dip(// = i) = E ^^'-/i(^-i+M = /iH^) = -*/f ■ 



(39) 



(40) 
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Now the matrix element ([37]) obtains a form: 

= (41) 

V 3 Wfi 

Comparing selection rules ([5^ in the dipole approximation with the selection rules ([77|) in the multipole approach, 
one can see that they do not impose any restrictions on the emission of photons. From eq. i3^) one can see that 
the wave function of photon in the dipole approximation has no any information about orientation of the vector k 
concerning vector r in any selected frame system. On the basis of such a fact it is logically to consider the emission 
of photons in the dipole approximation as isotropic. 



2. Differential matrix element 



Let US define the following differential components of the angular integral /dip by solid angle SI and by angle 6: 

-ln~ - sinedod^-ii ^ 2./i,-r,^„,(n.)ri,-,(n,), 

u— ±1 m 

rr (42) 

d hip {If) 



sinOdO , „ 

/J, — ±1 m 



But, in their calculations we shall assume that the selection rules ([59)1 are fulfilled. We find 



2 , 



dfl V 3 ^ ' Stt ' ' ' ' Stt 

2^ (43) 

didipilf) f dldip{lf) ■ v6 . 2/1 

— ' dip = —I — sm 9. 



sinOde J dn 4 



On their basis we shall define the following differential components of the matrix element: 

rfPdip ('/ = !) _ '/dip(l) dldip{lf) _ . Jdip{l) 



dfl Wfi dfl Stt Wfi 



sin^ I 



dpdipilf = l) _ Jdip(l) didipilf) _. Jdipi^) . 



I — '-^-^ ■ siv? I 



sinO dO Wfi sin 9 d9 4 Wfi 



(44) 



One can see that integration of such functions over the solid angle fl or over the angle 9 (with needed limits) gives the 
integrated matrix element pdip exactly. Comparing eqs. (j44p and (j4ip . we obtain connection between the differential 
and integral components of the matrix element in the dipole approximation: 

'^Pdip „ dpdip ^ • 2 /) fAl^\ 

~JFr = ■ ■ aM = Pdip ■ — sm^ 9. (45) 
dil sin9d9 Stt 

The differential matrix clement characterizes the angular distribution of emitted a-particles, which becomes anisotropic 
in result of emission of photons. 



G. Angular probability of emission of photon with impulse k and polarization e*^"^ 

I define the probability of transition of the system for time unit from the initial j-state into the final /-states, being 
in the given interval dvf, with emission of photon with possible impulses inside the given interval dvpf^, so (see ref.jSOjj, 
(42,5) § 42, p. 189; ref.;32'3, § 44, p. 191): 

dW = ^-^^^ ■dv = 2Tr \Ff,\^ 5{wf - w, + w) ■ dv, dv = dvj ■ dvph, dvph = t^-^t = ^ dwd^ph ^ ^^g^ 

where dvph and dvf are intervals defined for photon and particle in the final /-state, d^lph — d cos 9ph — 
sin9ph d9ph dipph, kph = w/c. Ffi is integral over space with possible summation by some quantum numbers of 
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the system in the final /-state. Such procedure is averaging by these characteristics and Ffi is independent on them. 
Then, interval dvf has only new characteristics and quantum numbers, by which integration and summation in Ffi 
was not fulfilled. Integrating eq. (|46l) over dw and substituting eq. ([2|) for Ffi, we find: 



dW 



m? 27r 



p{h,kf) 



dflph dvf, Wfi = Wi - Wf = f,^^ ■ ^^'^^ 



This is the probability of the photon emission with impulse k (and with averaging by polarization e^")) where the 
integration over angles of the particle motion after the photon emission has already fulfilled. 

To take direction of motion (or tunneling) of the particle after emission into account, I define the probability 
so: the angular probability concerning angle 6 is such a function, definite integral of which by the angle 9 with limits 
from to TT corresponds exactly to the total probability of photon emission |^7[ ): 

_«vM^ . £|4!!^L(,..,,) ^>--fe.%.«/) ,,,J. ,48, 

dilphd cos Of ZTTc^ TO^ 1^ a COS Of ) 

This probability is inversely proportional to normalized volume V. With a purpose to have the probability independent 
on V, I divide eq. (|48|) on flux j of outgoing a-particles, which is inversely proportional to this volume V also. Using 
quantum field theory approach (where v(p) = |p|/Po at c = 1, see [11], § 21.4, p. 174): 

j = n, ^(p,), V, = |v,| = ^ = (49) 

where Ui is average number of particles in time unit before photon emission (we have Ui = 1 for the normalized wave 
function in the initial i-state), v{pi) is module of velocity of outgoing particle in the frame system where colliding 
center is not moved, I obtain the differential absolute probability (while let's name dW as the relative probability): 

dP{^f,9f) _ dW{^f,ef) _ ^'//^' /p(fc^^fc^)^Zl^llMM+h.c.l. (50) 



dflphd cos 9 f dQphd cos 9 f hc^ki 2Tr c^ ki \ ' dcos9f 



H. Multipolar approach 



Let us find the bremsstrahlung probability in the multiple approach at the first values If = 1 and Iph = 1 . Starting 
from eqs. (jlSp and (|16p . and using the found differential and integral electrical and magnetic components (|30p and 
([3T|l . I calculate: 



Pi {h,kf) 
dpi {ki, kf) 



/^•{^(l,0)-^^/2.J(l,2)}, 



^ • |3J(l,l)-cos6i- V2J(1,0)- J(l,2)- (l-3f 



sin 9 d9 

and from eq. (|50p I obtain the absolute angular probability (i3q^ = 9f): 



■ sm 



d flph dcOS 'da-y 



^eff Wf^ 
8 TT nn? ki 



J(1,0)-^V2. J(l,2) 



J*(1,0) + ^J*(l,2) • (l-3sin2^?a' 
v2 ^ 



V2 



J*(l, 1) • cos 



(51) 



(52) 



+ h.e. > • sin -do 



I. Dipole approach 

Let us find the differential absolute probability ([50]) in dependence on angle 9 and the integrated absolute probability 
in the dipole approximation. Taking into account: 

f, ■ /2" Jdip(l) dpAip{lf = l) .\/6Jdip(l) . 2n 

Pdip {If = 1 = « — . ^—S^ = * —A sin ^' 

V 3 Wfi sin9d9 4 Wfi 
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I obtain: 



d flph d cos / 



7? p2 



& m? ki w fi 
E, 



^cff 



•iWfi 



Jdip (1) 



Jdip (1) 



sin^ 6*, 



(53) 



One can write: 



d ilph d cos 9 f 



37r-Pdip(0/) -sin^^. 



(54) 



Comparing angular dependence in such a result with the angular correlation obtained in the dipole approximation 
in (see eq. (9) in the cited paper) one can find that they coincide. 



J. Spectroscopic factor 



In order to take into account a non-unit possibility of formation of the a-particle in the state, from which it 
is further emitted outside, let us come to other problems of nuclear decays where this question has already been 
resolved. For example, in a problem of decay of nucleus through emission of proton we should take into account that 
the state which is occupied by the proton before its emission, is empty for the daughter nucleus. In order to obtain 
proper value for half-life we should divide it, calculated before by semiclassical approach (or others) directly, on the 
spectroscopic factor Sp^^ (see [1^). The spectroscopic factor can easily be calculated in the independent quasiparticle 
approximation (BCS), in which one assume that the ground state of odd-Z nucleus is one-quasiparticle state, while 
that of odd-odd system is two-quasiparticle configuration. In the BSC theory, the spectroscopic factor is given by 
5^*'^^ = u^, where is the probability that the spherical orbital corresponding to proton emitted is empty in the 
daughter nucleus [H, [s^l . For different proton emitters inclusion of the spectroscopic factors into formulas of half- lives 
improves agreement between calculated half-lives and their experimental values. Similar approach seems to be in the 
problem of a-decay, where the spectroscopic factor is successfully applied. 

On such a basis, let us modify the formulas for probability. For the multipole and dipole approaches from eqs. 
and (|53p we obtain the following modified absolute angular probabilities: 



dP. 



El+Ml 
jlt, 1 



dVlph dcosdn 



7,^ p2 



Wfi Ei 



(th) 



J(l,0)-^x/2. J(l,2) 



J*(1,0) 4- -^J*(l,2) • (l - Ssin^iJ^^j _ _L j*(i^i) . cosi9c«7 +h.e.| • sin^,?^^, 



(55) 



d riph d cos 9 
PdipiOf) 



E,, 



ryl 2 



Ei 



Jdip (i; 



(56) 



Jdip (1) 



One can see that inclusion of the spectroscopic factor Sa'^^ into formulas of the bremsstrahlung probability raises the 
spectra obtained before without it by eqs. (j52p and (|53p . However, in this paper we shall restrict ourselves only by 
the first preliminary calculations of the bremsstrahlung spectra for the deformed ^^^Ra nucleus. 



III. CALCULATIONS AND ANALYSIS 



In order to estimate efficiency of the model and accuracy, which it gives in determination of the angular absolute 
probability of the photon emission, I have calculated the spectra for the ^^"Po, ^^''Po, ^^^Ra and ^^''Cm nuclei. Here, 
the bremsstrahlung probability is calculated by eq. (15^ . The nucleus-a-particle potential is defined by eqs. (11)-(15) 
with parameters — by eqs. (16)~(22) in [23]. Q^-value is 5.439 MeV for ^lOpo^ 7.865 MeV for ^^^Po, 4.904 MeV 
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for ^^^Ra, 5.940 MeV for •^^'*Cm, according to ref. [s^ (see p. 63). The angle "da-^ between the directions of the 
Qf-particle motion (with possible tunneling) and the photon emission for ^^''Po and ^^^Ra is 90°, that is explained by 
configuration in experiments [l^, [l^Isl]. For ^^"Po I use the same angle but results could be easily generalized for 
other its values. For ^^''Cm I study two cases at 90° and 25° in order to include experiments [28*1 into the total picture. 
The first step is calculations of wave functions which should be obtained with high accuracy in order to obtain the 
convergent spectra. The wave functions of the a-decaying system in the states before and after emission of photon 
are presented in Fig. [1] (presentation is for the ^i°Po nucleus, energy of photon is 500 keV). One can see that in the 
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FIG. 1: Wave functions of the a-decaying system for the ^^"Po nucleus: (a) Imaginary part of the wave function in the initial 
i-state; (b) Real part of the wave function in the initial i-state; (c) The real wave function in the final /-state (after emission of 
photon); (d) Errors of the wave function in the final /-state appeared in current calculations of the spectra (for presentation, 
module of such a wave function is shown). 



initial state the wave function is complex. This provides fiux to be constant (in the internal region, in the region of 
tunneling and in the external region) and be directed outside. In current calculations of the spectra the wave function 
in the final state is real that allows to miss discontinuity at r = 0. In particular, in these figures (a, b, c) one can 
see behavior of these wave functions at boundary between tunneling and external regions. The last figure (d) shows 
errors which calculations give us, and one hope this allows us to obtain convergent pictures of the bremsstrahlung 
probability for all studied nuclei. 

According to analysis, the multipolar approach provides more accurate calculations of the bremsstrahlung spectra 
both in absolute scale (i. e. without any normalization on existed experimental data) and in a case when such 
normalization is used, in comparison with the dipole approach (see Appendix iDl for details). For example, let us look 
on Fig. [2] [left panel], where the spectra for the ^^°Po nucleus obtained by the multipole and dipole approaches in 
absolute scale are presented. In order to feel sensitivity of the spectra on change of the internal well of the a-nucleus 
potential (i. e. its shape in the internal region from r — > up to internal turning point), the following calculations 
are included into this figure: (1) the spectra without emission from such internal well (for the dipole approach this 
corresponds to a case of rectangular well in the internal region) , and (2) the complete spectra obtained concerning full 
a-nucleus potential with realistic nuclear component. It is clearly seen that for the dipole approach the influence of 
the shape of the potential in the internal region on the spectrum is stronger essentially while the multipole approach 
seems to be less sensitive to it. Such a peculiarity could be explained by more accurate use of far asymptotic region 
by the multipole approach, while the dipole approach reduces such a region. This leads to higher convergence of 
numerical integration over r in calculations of the spectra in the dipole approach. However, more accurate use of 
the far asymptotic region by the multipole approach provides essentially better agreement between obtained spectra and 
experimental data, while the dipole approach seems to do not calculate correctly here (at its present development). This 
important point adds the power of predictions in absolute scale to the multipolar model, in contrast to the dipole one. 
So, we shall use the multipolar approach for further estimation of the spectra for interesting nuclei and for analysis. 
In order to complete such a picture, I add calculations of the spectra for ^^"Po in the dipole approach for different 
values of the nuclear strength, presented in next Fig. O [right panel] . From here it is clearly seen that the spectrum is 
higher, if the nuclear strength is larger (i.e. the internal well of the a-nucleus potential is deeper). One can suppose 
that for multipole approach such dependence should be essentially smaller. 



A. Bremsstrahlung spectra for ^^''Po, ^'^^Ra and ^'^*Cm: comparison theory and experiments 

The best result in agreement between theory and experiment I have obtained for the ^^''Po nucleus (see the left 
panel in Fig. [31 here there is no any normalization of the calculated curve relatively experimental data). From figure 
one can see that for this nucleus the calculated spectrum by the proposed approach is in enough good agreement with 
the experimental data [l^ inside the region from 100 keV up to 750 keV. The calculated absolute probabilities of the 
bremsstrahlung in a-decay of the ^^^Ra nucleus and experimental data in [25j for this nucleus are presented in the 
central panel in Fig. [3] For this nucleus at low energies of the photons emitted the calculated spectra are located 
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FIG. 2: Bremsstrahlung probabilities in the a-decay of the ^'^"Po nucleus calculated in the absolute scale (solid line, red, is for 
the spectrum calculated by the multipole model; short-dash line, brown, for the spectrum calculated by the multipole model 
without emission from the internal region up to the internal turning point; dash line, green, for the full spectrum calculated by 
the dipole model; dash-dot-dot line, violet, for the spectrum calculated by the dipole model without taking the internal region 
up to the internal turning point into account; dash-dot line, blue, for the normalized spectrum calculated by the approach in 
[l| with normalization used in that paper): [left panel] difference between two spectra calculated in the dipole approach with 
inclusion of emission from the internal well and without such emission is not small, in contrast to the multipole approach; [right 
panel] a role of the nuclear strength Vq is clearly shown in calculations of the spectra in the dipole approach: one can see that 
the strength is larger, the spectrum is higher (for each spectrum the corresponding strength is defined as Vo = Vo ■ where Vb 
is defined by eq. (14) in ^31 a-nd additional new factor V has values 0,1,2,3,4,5 as shown in the figure). 



below experimental data, but for energies from 350 keV and higher I have obtained good agreement between theory 
and experiment. I add the calculated absolute probabilities for the ^"^^Cm nucleus, comparing them with the high 
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FIG. 3: The bremsstrahlung probability in the a-decay of the ^'^^Po, ^^'^Ra and ^^''Cm nuclei: red solid line is for the absolute 
probability calculated by the multipole model, dash-dot blue line for the normalized spectrum calculated by approach [25l ]. 
dash green line for the absolute probability calculated by formula (|55|) of the multipole model with taking the spectroscopic 
factor Sa^^ into account for the deformed nucleus ^^®Ra (we use approximated value Sa^^ = 0.2 from table I in [S^), scatter 
line for experimental data (data [H for ^"Po, data [H for ^^^Ra), data for ^''''Cm)). 

limit of errors of experimental data 0, (see the right panel in Fig. [3]). I calculate the bremsstrahlung spectrum 
for i!)aj = 25° of the angle 'da-y between direction of motion of the a-particle and photon emission which correspond 
to experiment [11], and I add another spectrum obtained for 'da-y = 90°. From this figure I see that both calculated 
curves are located enough close to the high limit of error of experimental data, and perhaps one can conclude that the 
agreement between theory and experiment is not so bad. This figure clearly demonstrates that the difference between 
higher and lower spectra could be explained by different values of this angle! Such explanation of two different spectra 
on the basis of one model is given for the first time, and one can suppose that by such a way a question discussed in 
[21I [23 ] concerning to the ^^°Po nucleus is explained also. 
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Q-decay dat 


a 


Bremsstrahlung probability, 1 / keV / decay 




Qa, MeV 


% 


-^1/2, a' ^'^'^ 


100 keV 


200 keV 


300 keV 


400 keV 


500 keV 


212 


7.987 


100.0 


3.0 E-2 


3.0 E-9 


8.1E-10 


2.7E-10 


9.5E-11 


3.5E-11 


218 


9.881 


100.0 


1.1 E-7 


7.5 E-9 


2.5 E-9 


1.0 E-9 


4.7E-10 


2.2E-10 


222 


8.164 


100.0 


2.8 E-3 


5.2 E-9 


1.3 E-9 


4.6E-10 


1.7E-10 


7.0 E- 11 


226 


6.487 


75.5 


2.5 E+3 


2.9 E-9 


5.6E-10 


1.3E-10 


3.5E-11 


9.4E-12 


228 


5.555 


72.7 


8.3 E+7 


1.8 E-9 


2.8E-10 


4.9E-11 


l.OE-11 


1.9E-12 



TABLE L Estimated values of the bremsstrahlung probability during a-decay of the Th nucleus and its isotopes 



B. Bremsstrahlung dependence on Qa and predictions of the bremsstrahlung probability during a-decay of 

isotopes of Th 

In Ref. [1^ it was reported about current investigations of bremsstrahlung accompanying the a-decay of the ^^*Th 
nucleus. Let us estimate the absolute bremsstrahlung probabihty for this nucleus on the basis of the proposed model. 
Results of such calculations are presented in Fig. [4l In calculations I use [11]: the angle between the directions of 
the a-particle motion (with possible tunneling) and the photon emission is 90°, Q^-value is 5.555 MeV. In Ref. (25j 
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FIG. 4: The predicted absolute bremsstrahlung probabilities in a-decay of the Th nucleus and its isotopes 



we explained the difference between the photon emission probabilities (both experimental and theoretical results) 
in the a-decay of ^^^Ra and ^^''Po (at first, dependence of the bremsstrahlung probability on the a-particle energy 
was analyzed in Ref. [2^): "The difference between the two sets of data can be attributed to the different structure 
of the two nuclei, which affects the motion of the a-particle inside the barrier. The ratio between the two sets of 
data of the photon emission probability dP/dE^ is strongly characterized by the different a-decay energy for ^^^Po 
(Ea — T.1 MeV) and ^^^Ra (Ea—4--8 MeV) concerning the shapes of the alpha-nucleus barriers for these nuclei." The 
difference between the a-particle energies for the decaying ^^"^Po and ^^^Ra nuclei is directly connected with different 
tunneling regions for these nuclei, which is directly connected with different contributions of the photons emission 
from tunneling and external regions, interference terms into the total spectra. And we obtained the property: The 
tunneling region is larger, the bremsstrahlung spectrum is smaller. The smaller values of the calculated total emission 
probability for ^^^Ra than the one for ^^^Po can be explained by a consequence of the fact that outside the barrier 
the Coulomb field (and its derivative respect to r) that acts on the a-particle in the case of ^^^Ra is smaller than in 
the case of ^^''Po because the external wide region results for the ^^''Po nucleus larger than for ^^^Ra and therefore 
the 7-emission probability for the ^^''Po nucleus is bigger. In Fig. 2] one can see the demonstration of this property 
for isotopes of Th. In Table [J it is shown how the bremsstrahlung probability depends on Q^-value of the nucleus for 
different energies of the photons emitted. 
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C. How much is the bremsstrahlung probabiHty changed in dependence on numbers of protons and 

neutrons of the a-decaying nucleus? 

Let us analyze how much the probabiUty of photons emitted could change if to change number of protons or neutrons 
of the nucleus which decays. In order to made such analysis clearer, I have fixed Q-value and have calculated the 





of Pf° 



FIG. 5; Distribution of the absolute bremsstrahlung probability in the a-decay on numbers of protons and nucleons of the 
decaying nucleus (used data: Qa = 5.439 MeV for ^^°Po): left panel is for the probability of the photons emitted with energy 
E~i = 100 keV, right panel for the probability of the photons emitted with energy = 300 keV 

probabilities in some region of numbers of protons and nucleons (it is convenient to use angle between directions of 
the emission of photon and motion of the a-particle equal to 90° , as results are supposed to be similar for any other 
angles) . In the next Fig. [S] a distribution of the absolute probability in dependence on the numbers of protons and 
neutrons of the a-decaying nucleus is presented for 100 keV, 300 keV and 500 keV of the emitted photons. One can see 
that at increasing of the number of neutrons the probability is changed a little and monotonously while at increasing 
of the number of protons the probability decreases stronger essentially and monotonously (such proton dependence 
could be a quantitative demonstration of results of 



D. Angular emission of photons 



Now let us analyze how much the bremsstrahlung probability is changed in dependence on the angle between 
directions of the motion of the particle and emission of photon. For analysis we consider the ^^"Po nucleus. In Fig. [5] 
the bremsstrahlung probabilities calculated by the multipole model at different values of this angle are presented. In 
particular, one can see that spectrum for the angle 25° is located enough far below then the spectrum for the 90°. Now 
these calculations confirm (at first time) arguments proposed for explanation of difference between experimental data 
[28l | and [11] obtained for these angles, correspondingly, and they have resolved discussions in [HI, [l!]. In the next 
Fig. [7] distribution of the bremsstrahlung probability on the angle 'da-y between directions of the motion of the particle 
and emission of photon for the ^^"^Po nucleus is presented. In particular, one can see that the angular distribution 
of the probability is formed mainly by the first integral J(l, 0) while next two integrals J(l, 1) and J(l, 2) give very 
small contribution to the total angular spectrum. 
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FIG. 6; The absolute bremsstrahlung probability in the a-decay of Po in dependence on the angle 19^7 between directions 
of the motion of the a-particle and the emission of photon calculated by the multipole model 




FIG. 7: Distribution of the bremsstrahlung probability on the angle between directions of the motion of the particle and emission 
of photon for the ^'^"Po nucleus: (a) angular probability at 300 keV, 350 keV, 400 keV, 500 keV; (b) angular probability at 
300 keV in large scale: here one can see difference between two spectra. 



E. Formula of the bremsstrahlung probability in the a-decay 

Let us restrict ourselves by only one nucleus and try to write such formula for it. After preliminary estimations of 
the spectra for different nuclei, I propose the following form (further, let's use only -da-^ — 90°); 

, f dPparam(w; Oq... 04, ni... 724)1 f ^Iff^A , „, , a2 03 ^4 ^r^N 

In { — —rr: ; 7, } = In i r — -rr, — > + flo - ai w ^ H , (57) 

where ao . . . 04 and ni ... 7^4 are unknown constants which do not depend on energy of the photon emitted and are 
changed for the different nuclei. These constants reflect "structure" of the a-decay for the studied nucleus. Therefore, 
they should depend on Qa, ZeS, Zd, Ad of this nucleus. 

With a purpose to find parameters oq . . . a4 ni . . .n^ for the selected nucleus, I shall introduce the following char- 
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a-decay data 


parameters 


No. 


Ad 


AT 


Zd 




Qc MeV 


(mm) 
% 


Iparamj 

"o 




IparamJ 
"l 


1 


224 


6.073177 


88 


0.42105 


5.555 


10.2 


10.20083 


0.0154 


0.01531749 


2 


224 


6.073177 


88 


0.42105 


10.0 


11.2 


11.20020 


0.0069 


0.00681732 


3 


102 


4.672328 


50 


0.03774 


10.0 


6.3 


6.30084 


0.00475 


0.00440210 


4 


262 


6.398827 


107 


0.36090 


10.0 


10.9 


10.90000 


0.008 


0.00799993 



TABLE II: Parameters ao and ai for ^^^Th, 
calculated by the method of minimization, c 



'Te and nucleus with Ar, 



266 and Zt, 



(par am) 
'0 



and a I 



(param) 



109 (qq™'"' and a^™'"' are parameters 



are parameters calculated by formula ((6TJ) 



acteristic: 



<7 {tti, Hi) 



\ 



AP {w; Ui, n.i) = In 



^max ^min 



dPaiodcl{w) 



d flph d cos 6 



AP {w; ai, Hi) dw, 

dPparaiYiiw] ao . . . a4, ni . . . 714) 



(58) 



In 



dVlph d cos Of 



where dPmodoi and dPparam are the bremsstrahlung probabiUties calculated by the multipole model and by the formula 
([57)) . correspondingly. The a at selected set of parameters is smaller, the curve dPparam obtained by formula (|57)l is 
closer to the spectrum dPmodci calculated by the multipole model. I. e. the best description of the bremsstrahlung 
spectrum for the studied nucleus by formula (|57p should be obtained at such choice of the parameters oq . . . 04 and 
ni ... 714 where a is minimal. For convenience, I call this approach for determination of parameters for the selected 
nucleus as method of minimization. So, using the method of minimization, for the ^^*Th nucleus I obtain the following 
values (wmin = 50 keV and w^ax = 900 keV are used): 



ni — 1, 7i2 = 0.5, rt3 = 1, ?i4 = 2, 
ao = 10.8, ai = 0.007, as = 10, ag = 10, a4 = 1. 



(59) 

The curve calculated by formula ([57]) at choice ([55]) of parameters turns out to be located extremely close to the 
bremsstrahlung spectrum for '^^^Th! From here I conclude that the bremsstrahlung probability for arbitrary one nucleus 
can be approximated by formula ^57\ l with very high accuracy inside the energy region up to 1 Me V. Estimations of 
parameters for other nuclei show that it is possible to describe the bremsstrahlung spectra with enough high accuracy 
for different nuclei using different values of the oq and ai parameters only, while the ni ... 77.4 parameters and even 
the 02, 03, 04 parameters are fixed. 

In this paper, I shall define the tii . . .774, a2, as and a4 parameters for the different nuclei by (|59p and try to find 
out how ao and ai can be described. Assuming dependence of oq and ai on Q, Ad and Zd to be linear, I propose the 
following formula: 



flo (Q, Ad, Zd) = boo + 5oi Q + feo2 Ad + 603 Zd, 
ai{Q,Ad, Zd) = bio + biiQ + bi2Ad + bi3Zd, 



(60) 



where new unknown parameters feoi and bu {i = 0, 1, 2, 3) have been introduced which do not already depend on Q, 
Zd and Ad- The simplest way is to find 601 and 61 1. For such calculations we needs in two nuclei with equal Zd and 
Ad values and different Q-values. Let's consider the ^■^^Th nucleus. For it I calculate the bremsstrahlung spectrum 
on the basis of the multipole model at two different Q-values (I use: Qi = 5.555 MeV and Q2 — 10 MeV), and then 
I obtain the ao and oi parameters using the method of minimization above. Results are presented in Table [TTl in the 
first two strings with numbers 1 and 2. 

In order to find next unknown parameters bij (i — 1, 2, j — 0, 2, 3), it needs to consider nuclei with the different Ad, 
Zd numbers at the same Q- value. To achieve accuracy as high as possible, we shall use the previous nucleus and two 
other nuclei with the largest difference between Ad (and between Zd). Let's use Table in (S^, from here we select: 
^"^Te and nucleus with Ap = 266, Zp = 109. I calculate the bremsstrahlung spectra at Qa-value equals to 10 MeV 
using the multipole model and then I find ao and ai for them using the minimization method. Results are presented 
in Table HTl in the next two strings with numbers 3 and 4. Using data of the table HIl I calculate unknown boi and bu, 
then I find the final form of oq and ai on Q, Ad and Zd {Q is used in MeV): 



aQ{Q,Ad,Zd) = 4.60202 + 0.22497- g + 0.11956- Ad - 
ax {Q,Ad,Zd) = 0.0204108- 0.0019123- g + 1.086956 



0.25492 - Zd, 

10-6 - Aa + 6.0068649 ■ 



10- 



Za 



(61) 
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and the bremsstrahlung formula ([57]) has transformed into such: 

^^f dPf'+''\w,9f = 90°) \ _ f ^eff^^ \ , I 10 I 10 I 1 (62) 

\ d flph d cos 9f J \ 8 TT TO^ /Ci J -y/w w w"^ 

For four studied nuclei I have obtained very small difference between the probability calculated by the model above 
and the curve calculated by formula ([5^ with parameters for the energy of photons emitted up to 1 MeV (less 
then 1 percent). /. e. we have described the bremsstrahlung spectra inside the energy region up to 1 MeV for four 
different nuclei (with such long maximal distance between their numbers A^) with very good accuracy by only one 
this formula with parameters calculated only on the basis of values Ad, Z^, Qa! It turns out that description of the 
bremsstrahlung spectra for all nuclei inside region Ad — 107.. 262 at different Zd by the formula ((62)) with parameters 
(|6ip is not such accurate but enough good also (see curves in Fig. [5] obtained by formula). However, one can improve 
further such approximation essentially for the "problem" nuclei, if to pass from the linear dependence (|60p of the 
bremsstrahlung probability on the Ad and Zd values to harmonic one. 
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FIG. 8: The bremsstrahlung probabUity in the ct-decay of the ^^'^Po, ^^'^Po and ^^^Ra nuclei (i?o 



90°): the absolute 



probability calculated by the multipole model (red solid line), experimental data (scatter line, data [2J] for ^^^Po, data ^ for 
^^"Po and data ^25] for ^^^Ra) and curve calculated by formula (|57p with (|6ip (dash blue line). 



IV. CONCLUSIONS 



In the paper the model of the bremsstrahlung accompanying the a-decay is presented where emphasis is given on 
development of unified angular formalism of the dipole and multipole approaches. Effectiveness of the model and 
accuracy of calculations of the bremsstrahlung spectra are analyzed in their comparison with experimental data for 
the 2iopo, 2i4po, 226j^j^ ^^^^ 244^^^^ jv^^^g ^j^g following. 

• The multipolar model is the most motivated from the physical point of view, it is the richest in obtaining useful 
information about angular emission of photons during a-decay, their results obtained for the ^^°Po, ^^^Po and 
^^^Ra nuclei both in the absolute scale and with normalization on experimental data [1, HJ, [2^ are in the best 
agreement with these experimental data in comparison with any other known models and approaches. 

— The best result has been obtained in agreement between the calculated absolute probability of the 
bremsstrahlung emission for the ^^*Po nucleus and the experimental datajl^l for this nucleus inside the 
region of photons energies from 100 keV up to 750 keV (see Fig. [3] (a), Qa = 7.865 MeV, the angle daj 
between the directions of the a-particle motion and the photon emission is 90°). 

— The calculated absolute probabilities of the bremsstrahlung emission in a-decay of the ^^°Po and ^^^Ra 
nuclei for low energies of the photons emitted are located below experimental data ^ and [255, lo^' 
energies from 350 keV and higher I have obtained good agreement between the model and experiment (see 
Fig. [3] (b) and (c), Q„ = 5.439 MeV for ^lOpo and = 4.904 MeV for ^^^Ra, i}a-y ^ 90°). 

— The spectrum for the ^''^Cm nucleus obtained at -da-y — 25° is found to be in satisfactory agreement with 
high limit of errors of experimental data ^28,] . 

• In frameworks of the dipole approach, the emission of photons from the internal region before the barrier is not 
small and gives visible contribution into the total spectrum. However, the multipolar approach seems to be less 
sensitive to such emission. Such a peculiarity could be explained by more accurate use of the far asymptotic 
region by the multipole approach, while the dipole approach reduces this region and, so, has higher convergence 
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in numerical integration of the spectra over r. However, more accurate use of the far asymptotic region by 
the multipole approach provides essentially better agreement between obtained spectra and experimental data, 
while the dipole approach seems to be fail here (at its present development). This important point adds the 
power of predictions in absolute scale to the multipolar model, in contrast to the dipole one. This demolishes 
published progress in agreement between experimental data and the spectra calculated in the dipole approach 
(where such emission from the internal region was neglected) if it was affirmed as obtained in the absolute scale. 

• A hypothesis about explanation of difference between two experiments [s^ and on the basis of different values 
25° and 90° of the angle between direction of motion of the a-particle and emission of photons proposed in 
discussions [ll], [S^l, has been confirmed (at the first time). 

• The unified formula of the bremsstrahlung probability (in the absolute scale) during the a-decay of the arbitrary 
nucleus, which is directly expressed through the Qa-value and numbers Ad, Zd of protons and neutrons of this 
nucleus, has been constructed. Inside region of the a-active nuclei from ^'^^Te up to the nucleus with numbers 
of nucleons and protons Ap — 266 and Zp — 109 (this region is taken from [35|) with energy of the photons 
emitted from 50 keV up to 900 keV satisfactory agreement has been achieved between the spectra, obtained 
on the basis of the multipole model (where duration of calculations for one selected nucleus is up to 1 day), 
and the bremsstrahlung spectra obtained on the basis of the proposed formula (where duration of calculations 
is about some seconds using the same computer). However, analyzing results for the ^^'^Po, ^^°Po and ^^^Ra 
nuclei, this formula is found to give more accurate spectra relatively experimental data (up to 500 keV), then 
spectra obtained in the dipole approach without taking the emission from the internal region before the barrier 
into account. 
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APPENDIX A: CLEBSCH-GORDAN COEFFICIENTS 



I define Clebsch-Gordan coefficients, according to Table I in ref. [3l[ (see p. 317) and find: 



(Olll 2,-1, 1) = 0, (ml 2,-1, 1) = 0, (21l|2,-l, 1) 



5' 



(011|0, 1, 1) = (111|0, 1, 1) = -Ji, (211|0, 1, 1) = J^, 



(011|0,-1,-1) 



(lll|0,-l,-l) 



\, (211|0,-1,-1) = J1, 



(Al) 



(Olll -2, 1,-1)^0, (ml - 2, 1,-1) = 0, (2111-2,1,-1) 



APPENDIX B: COEFFICIENTS C] 



We define the coefficient C,'"C „ so: 



l{2lf + l){2n+l) (If-iy. (n-|m + ^'|)! 



X 



327r (^/ + 1)! (n+|m + ^'|)! 



At Z/ = 1, Iph = 1 and n = we liave: 

m = — /i' = ±1 

and tfie coefficient CT^i" „ is: 



At If — 1, Iph ^ I and n — 1 the property (|B2[) is fulfilled and we obtain: 



• (llll 0,^',-m). 



647r 

At Z/ = 1, Zp/i = 1 and n ^ 2 the property (|B2[) is not fulfilled and 



^"V _ ^ / 15 (2- |m + A^^|)! , 



Rewrite these coefficients at different m — ±1 and n' — ±1: 



-1-1 _ 1. 5 . .^ ^-11 1 /15 



1 [15 . 1 



^iV = -gVv^'^^l^-^'-^)' - 16 -('''I 

Substituting here values (|A1|1 for the Clcbsh-Gordan coefficients, we find: 

-1-1 _ n r^-ii _ 1 . / ^ /^i-i _ 1 . / ^ ^11 



<^iio -0, Clio --8'V2^' '^"•^ ""8 V 2^' ^"0-°' 

8 V 27r' 1" 8 ■ V 2^' 



/^-l-l _ n /-i-ll _ . . / /^l-l _ _ . , / fill _ n. 

•-^111 — "-"111 ~ o Vo_' ^111 ~o V ^111—"' 



^-1-1 _ 1 / 3 ^-11 _ -'^ / ^1-1 _ -'^ / ^ r<ii _ ^ 

" 16 V2^' ""8V2^' ^"2"~8V2^' ^"^"16 



APPENDIX C: FUNCTIONS f[^'^ (9) 

Let us consider the function {9): 



fZ'iO) = PA (cos 0) pi (cos 0) 4'"+'^' I (cos ( 
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At Z/ = 1 and n = 0, 1, 2 we obtain: 

Pi (cos 61) Pi (cos 61) P(l"+'^' I (cos 9) , 

Pi (cos 9) Pi (cos 9) pI"'^''' I (cos 9), (C2) 
pi(cos6i) pi(cos6i) pl"+'''l(cos6l). 

At different m = ±1 and /i' = ±1 we find: 



Jio y^i 


= Pi {cos 9 

— l^COS (7 

= Pi {cos 9 
= Pi {cos 9 


) Pi {cos 9] 

) Pi {cos 9] 
) Pi {cos 9] 


PS{cos9] 

pf{cos9] 
pf{cos9] 


= 0, 
— sm (7, 
= sin^e, 
= 0; 




fim 


= Pi {cos 9 
= Pi {cos 9 
= Pi {cos 9 
= Pi {cos 9 


) Pi {cos 9] 
) Pi {cos 9] 
) Pi {cos 9] 
) Pi {cos 9] 


P^{cos9] 
pf{cos9] 
pf{cos9] 
pf{cos9] 


= 0, 

— sin^ cos 9, 
= sin^ cos 9, 
= 0; 






= Pi {cos 9 
= Pi {cos 9 
= Pi {cos 9 
= Pi {cos 9 


) Pi {cos 9] 
) Pi {cos 9] 
) Pi {cos 9] 
) Pi {cos 9] 


Pi {cos 9] 
pf{cos9] 
pf{cos9] 
Pi {cos 9] 


= Ssin'^e, 
= i sin^ 6* (3 cos^ 6 
= i sin^ 9 (3 cos^ 6 
= 3sin^6'. 


-1), 
-1), 



(C3) 



APPENDIX D: DIPOLE APPROACH VERSUS MULTIPOLE ONE: COMPARATIVE ANALYSIS AND A 

ROLE OF THE INTERNAL WELL OF THE POTENTIAL 

After appearance of the fully quantum approach proposed by Papenbrock and Bertsch in [l^ where wave function 
of photons was used in the dipole approximation, further fully quantum approaches have been developed mainly on 
its basis. In particular, formula (/|p|i) = ih {f\drV\i) / Ej proposed in this paper for transformation of the matrix 
element of the photon emission and increasing essentially its convergence in calculations without visible decreasing 
of accuracy, becomes very popular. This is essential point which attracts many researchers to study this problem in 
fully quantum approach. Published results with agreement between spectra calculated by such dipole approach and 
experimental data look to be well. In particular, consequence "The present high precision data clearly demonstrate 
the failure of a classical Coulomb acceleration calculation (see, e.g., [9,14]) to describe the bremsstrahlung emission 
in a decay, and rules out theoretical suggestions put forward by authors of Refs. [6,7,9]" in ^ (see p. 4, reference 
in cited paper) has been considered further as a fact that the dipole approach is the most accurate in description of 
experimental data and perspective from all fully quantum approaches in further study of bremsstrahlung emission 
during a-decay. However, let us clarify how much such approach is accurate in description of experimental data in 
comparison with the multipole one. 

1. Spectra in absolute scale 

Usually, authors do not mention whether they calculate normalized or absolute probability in the dipole approach. 
However, the spectra calculated in the absolute scale are sometimes noted to have a main progress of such a way. 
But if to suppose that these results were obtained without any normalization on experimental data, then I meet the 
following problem. Up today, all published calculations in the dipole approach have been based on the a-nucleus 
potential where rectangular well inside the internal region before the barrier was used. According to ()33|) . we obtain 
directly null contribution of the emission from this internal region into the total spectrum. But such results would be 
reliable if the real emission from the internal region is very small. Only in such a case more accurate realistic shape of 
the well could be neglected. But I find that this is not so. Let us look on the left panel of Fig. [9] where the spectra for 
the ^^^Po nucleus calculated by the dipole and multipole approaches without any normalization on experimental data 
are presented. Here, one can see that difference between the full spectrum obtained by the dipole approach (see dash 
line, green, in this figure) and the spectrum obtained by the dipole approach formed in the spatial region starting from 
the internal turning point (see dash-dot-dot line, brown, in this figure) is not small (in calculations for ^^''Po such data 
are used: Qa = 5.439 MeV, internal turning point is rtp.int = 8.829 fm, external turning point is rtp,ext — 44.007 fm, 
tunneling length is Artun = 35.177 fm)! One can find that complete neglect of the emission from the internal region 



/iT>) - 
/iT>) = 
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E(keV) EJkeV) 

FIG. 9: Bremsstrahlung probability in the a-decay of the ^^°Po nucleus (solid line, red, is for the spectrum calculated by 
the multipole model; dash line, green, for the full spectrum calculated by the dipole model; dash-dot-dot line, brown, for the 
spectrum calculated by the dipole model without taking the internal region up to the internal turning point into account; 
dash-dot line, blue, for the normalized spectrum calculated by the approach in with normalization used in this paper): (left 
panel, a) Absolute probabilities: here one can see that difference between two spectra calculated in the dipole approach with 
inclusion of the internal region and without it is not small that confirms a real non-small influence of this nuclear region on 
the total spectrum: (right panel, b) Normalized probabilities: here one can see that normalized spectra in multipole approach 
and in approach [J] are found in better agreement with experimental data 0] in comparison with normalized results obtained 
in the dipole approach (in contrast with conclusions of [H). 

improves visibly agreement between the dipole calculations and experimental data. However, these two spectra are 
essentially farther to experimental data in comparison with the spectrum obtained by the multipole model without 
any normalization (see solid line, red, in this figure). Now this confirms a real importance of inclusion of the realistic 
shape of the well before the barrier into calculations of the spectra in the absolute scale. But if the spectra obtained 
previously in the dipole approach were in good or the best agreement with experimental data (authors are assumed 
to have catched a successful normalization in theory), then after taking the emission from the internal realistic well 
into account such spectra should be displaced below (to the distance similar to distance between two curves in the 
dipole approach in Fig. [5]). As a result, this demolishes all published progress in agreement between experimental data 
and calculations in the dipole approach. Now physical motivations of such progress in agreement between theory and 
experiment have became unclear if they were affirmed to be obtained in the absolute scale. It has became clear which 
difficulty authors of [lB| had, who gave the first predictions for the ^^°Po nucleus before the first experiments [36j . 

2. Spectra normalized on experimental data 

If to suppose that the spectra have obtained with normalization on experimental data, then we must compare them 
with results obtained in the multipolar approach with normalization also. If in eq. (j52p for the absolute probability 
in the multipolar approach to restrict ourselves by the first most important integral J(1,0) only in comparison with 
other two integrals J(l,l) and J(l,2), then we obtain (with possible normalization) the spectrum calculated by the 
model [l| exactly. Now exact coincidence has been established between the spectra obtained in the approach [l[ and 
the spectra obtained by the multipolar model, where a difference is explained by taking non-zero magnetic component 
at J(l, 1) into account in the total matrix element in cqs. (jSTjl (J(l,2) is smaller then J(1,0) by 1-2 orders from 50 
keV up to 1 MeV of the photons). The matrix element beyond the dipole approximation through expansion of the 
wave function of photons in spherical waves (at first, proposed in 23]) and realistic form of the a-nucleus potential 
realized in the model fl| give us a small difference between the spectrum calculated for the ^-'^^Po nucleus by such 
a way and the spectrum obtained for this nucleus in the dipole approach with formula (/|p|«) = ih {f\drV\i} / 
(see lines 6 and 7 in Fig. 1 in [l] and comments here). As latter experiment [Sj for ^^°Po showed, this predicted 
peculiarity increases a little coincidence with experimental data of such paper in comparison with our calculations 
in the dipole approach above and with presented results (at lower energies of photons) in the fully quantum dipole 
approach in this paper also. While authors of 0] did not comment such comparison, did not present any grounds 
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of their affirmation about groundlessness (ruling out) and accuracy of the approach the complete comparative 
consideration of these experimental data and all these spectra breaks down such affirmation (see Fig. 3 in ,.25.] and 
discussions). The same logics is applicable again for the multipole model directly. Let us look to the right panel in 
Fig. O Here, one can see that the spectrum calculated by the multipole model (see solid line, red, in this figure) looks 
to be a little more successful in description of the experimental data ^SJ in comparison with the spectrum calculated 
by the dipole approach (see dash line, green, in this figure). Here, our discussed result by the approach [25| is included 
also (see dash-dot line, blue, in this figure). Results in descriptions of the newest experimental data [2, lU, for 
the ^^''Po and ^-^^Ra nuclei in direction of the approach Jj look enough good also, where we have been achieving 
agreement between theory and experiment up to 765 keV (in |^ energy region of photons emitted is up to 500 keV). 

Irrespective of accuracy which the approach could give to researcher, more important point consists in physics 
which such approach has. While calculations of spectra in all other published approaches up today [gI. [sI. [qI. [TqI. l20l [27| 
are based on the potential of interaction between the a-particle and the daughter nucleus, which in the region of 
nuclear forces has a form of rectangular well and its depth is determined for each selected nucleus separately and 
absolutely does not take the real nuclear shape into account, our model and calculations starting from [l| have been 
based on the unified global realistic nuclcus-a-particle potential which parameters are defined by Qa-value of the 
a-decay, protons and neutrons numbers only for the studied nucleus. Such potential is constructed on the basis 
of analysis of 344 a-decaying nuclei in frameworks of one unified model UMADAC of the a-decay and a-capture 
[40I I4H (see also 37]) where errors with experimental data turn out to be the smallest in comparison with other 
known models this year. This unification of the potential combines naturally with our formalism for calculations of 
the absolute angular (and integral) probabilities of bremsstrahlung during the a-decay. This allows to calculate the 
absolute bremsstrahlung probability for arbitrary nucleus using Q^-value, proton and neutron numbers only as input 
data. After angular realization of the multipolar approach in calculation of the matrix elements, such a way allows 
to take a deformation of the a-decaying nucleus into account [2] , which turns out to be not small (that confirms a 
real importance of accurate determination of the a-nucleus potential in the region of the nuclear forces) and it could 
be extracted from the bremsstrahlung spectra (what other published approaches in the dipole approach are not able 
to study in current their stage) . This fact underlines perspective of further investigations of the multipole approach. 
One can hope it should allow to open new independent experimental ways to obtain new information about a-decay. 

Argument presented above contradict to the affirmation of authors in [8] about groundlessness of the approach [l| 
and results obtained by such approach. So, among all variety of existed models and approaches the multipolar model 
is the most motivated from the physical point of view, it is the richest in obtaining useful information about emission 
of photons during a-decay, their results are in the best agreement with experimental data existed. 



APPENDIX E: TRANSFORMATIONS OF THE MATRIX ELEMENT IN A GENERAL CASE 



The matrix element of emission p2|) in the dipole approximation can be transformed as [1^ 

dU{v) 



(/IpN) = 



ih(f 



dv 



(El) 



E,-E 



f 



In the resulting expression the integrant function obtains additional factor l/r at enough far r (in result of Coulomb 
type of the a-nucleus potential here), that in the asymptotic region increases convergence of numerical integration of 
the total matrix element over r. So, it could be useful to generalize such a transformation for the multipole approach. 
However, in a general case the calculation of the matrix element looks to be more complicated and the following 
theorem can be applied. 



3(a) gi (kr-u)t) 



Theorem. Let us consider transition of the a-decaying system from the initial i-state into the final f -state in result 

of emission of photon. If the vector potential of photons was used in such a form A (r) 
then the matrix element of such transition i ^ f would be written as 

dU{r) - ^2 



zhlf 



(/|A*p|z) = - 



dr 



m 



{' 



d_ 
dr 



d_ 
dr 



he 

2m 



(E2) 



where U (r) is the a-nucleus potential, m is reduced mass, E^ and Ef are energies of the a-decaying system in the 
initial i-state and in the final f -state, k is wave vector of photon emitted, k = |k| is its wave number. 
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Now taking into account 

dU(r) 

[H,p]=in^^, (ElO) 

we obtain ((E2)) . □ 

Physical sense of the formula (IE2p could be explained by the following. The first item in this formula seems to 
have higher convergence in the asymptotic region at numerical integration over r. It is supposed to give major 
contribution into the total matrix element while the second item allows us to analyze corrections after taking more 
accurate estimation of the emission of photons in the far asymptotic region into account (which is smaller usually). 



S. P. Maydanyuk, V. S. OIkhovsky, Europ. Phys. Journ. A28 (3), 283-294 (2006), nucl-th/0408022' 

S. P. Maydanyuk, V. S. OIkhovsky, G. Giardina, G. Fazio, G. Mandaglio, M. Manganaro, Nucl. Phys. A823 (1-4), 38-46 
(2009). 

M. I. Dyakonov, I. V. Gornyi, Phys. Rev. Lett. 76 (19), 3542-3545 (1996). 
M. I. Dyakonov, Phys. Rev. C60, 037602 (1999), nucl-th/9903016, 

C. A. Bertulani, D. T. de Paula and V. G. Zelevinsky Phys. Rev. C60 (3), 031602 (1999) (4 pages), nucl-ex/9812009", 
N. Takigawa, Y. Nozawa, K. Hagino, A. Ono and D. M. Brink, Phys. Rev. C59 (2), R593-R597 (1999), nucl-th/980900l] 
J. Kasagi, H. Yamazaki, N. Kasajima, T. Ohtsuki and H. Yuki, Bremsstrahlung emission in a-decay and tunneling motion 
of a-particle, Journal of Physics G: Nuclear and Particle Physics 23, 1451-1457 (1997). 

H. Boie, H. Scheit, U. D. Jentschura, F. Kock, M. Lauer, A. I. Milstein, I. S. Terekhov, and D. Schwalm, Phys. Rev. Lett. 
99, 022505 (2007), arXiv:0706.2109 

U. D. Jentschura, A. I. Milstein, L S. Terekhov, H. Boie, H. Scheit, and D. Schwalm, Phys. Rev. C77, 014611 (2008). 
S. Misicu, M. Rizea and W. Greiner, Journ. Phys. G 27, 993-1003 (2001). 
W. van Dijk and Y. Nogami, Few-body systems Supplement 14, 229-232 (2003). 
O. Serot, N. Carjan and D. Strottman, Nucl. Phys. A569, 562-574 (1994). 
W. van Dijk and Y. Nogami, Phys. Rev. Lett. 83, 2867-2871 (1999). 
W. van Dijk and Y. Nogami, Phys. Rev. C65, 024608 (2002) (14 pages). 



B. Ivlev and V. Gudkov, Phys. Rev. C69, 037602 (2004) (4 pages), nucl-th/0 307012 
M. Ya. Amusia, B. A. Zon, and I. Yu. Kretinin, JETP 105 (2), 343-346 (2007^^ 



V. V. Flambaum and V. G. Zelevinsky Phys. Rev. Lett. 83, 3108-3111 (1999), 'nucl-th/9812076" 
L S. Batkin, L V. Kopytin and T. A. Churakova Internal bremsstrahlung accompanying a decay, Yad. Phys. (Sov. Journal 
of nuclear physics) 44, Iss. 6 (12), 1454-1458 (1986). 

T. Papenbrock, G. F. Bertsch, Phys. Rev. Lett. 80 (19), 4141-4144 (1998), nucl-th/ 9801044 

E. V. Tkalya, Zh. Eksp. Teor. Fiz. 116, 390 (1999) [Translation: Sov. Phys. JETP 897T999)T08] ; E. V. Tkalya, Phys. Rev. 
C60, 054612 (1999). 

N. V. Eremin, G. Fazio and G. Giardina, Comment on "Bremsstrahlung in a-decay of^^'^Po: do a-particles emit photons 
m tunneling?" Phys. Rev. Lett. 85 (14), 3061 (2000). 

J. Kasagi, H. Yamazaki, N. Kasajima, T. Ohtsuki and H. Yuki, Replay on Comment on "Bremsstrahlung in a-decay of 
^^°Po: do a-particles emit photons in tunneling?" Phys. Rev. Lett. 85 (14), 3062 (2000). 

S. P. Maydanyuk, S. V. Belchikov, Prob. At. Sci. Tech.. Ser.: Nucl. Phys. Inv. (44) 5, 19-21 (2004), nucl-th/0404013[ 
G. Giardina, G. Fazio, G. Mandaglio, M. Manganaro, C. Sacca, N. V. Eremin, A. A. Paskhalov, D. A. Smirnov, S. P. May- 
danyuk, and V. S. OIkhovsky Europ. Phys. Journ. A36 (1), 31-36 (2008). 

G. Giardina, G. Fazio, G. Mandaglio, M. Manganaro, S. P. Maydanyuk, V. S. OIkhovsky, N. V. Eremin, A. A. Paskhalov, 

D. A. Smirnov and C. Sacca, Mod. Phys. Lett. A23 (31), 2651-2663 (2008), arXiv:0804.2640 

S. P. Maydanyuk, V. S. OIkhovsky, Prog. Theor. Phys. 109 (2), 203-211 (2003), n ucbth/04040 90l, 

W. So and Y. Kim, Journ. Korean Phys. Soc. 37 (3), 202-208 (2000). 

J. Kasagi, H. Yamazaki, N. Kasajima, T. Ohtsuki and H. Yuki, Bremsstrahlung in a-decay of^^'^Po: do a-particles emit 
photons m tunneling? Phys. Rev. Lett. 79 (3), 371-374 (1997). 

N. N. Bogoliubov and D. V. Shirkov, Kvantovie polya (Quantum fields theory) (Nauka, Mockva, 1980), p. 320 — [in 
Russian] . 

L. D. Landau and E. M. Lifshitz, Kvantovaya Mehanika, kurs Teoreticheskoi Fiziki (Quantum mechanics, course of 
Theoretical Physics), Vol. 3 (Nauka, Mockva, 1989) p. 768 — [in Russian; eng. variant; Oxford, Uk, Pergamon, 1982]. 
J. M. Eisenberg and W. Greiner, Mehanizmi vozbuzhdenia yadra. Electromagnitnoie i slaboie vzaimodeistviya (Excitation 
Mechanisms of Nucleus), Vol. 2 (Atomizdat, Moskva, 1973) p. 348 — [in Russian; Engl.: Excitation mechanisms of the 
nucleus. Electromagnetic and weak interactions (North-Holland publishing company, Amsterdam- London, 1970)]. 
V. B. Berestetsky, E. M. Lifshitz and L. P. Pitaevsky, Kvantovaya Elektrodmamika, kurs Teoreticheskoi Fiziki (Quantum 
electrodynamics, course of Theoretical Physics), Vol. 4 (Nauka, Mockva, 1989) p. 704 — [in Russian; eng. variant: Oxford, 
Uk, Pergamon, 1982, 652 p.]. 

S. Aberg, P. B. Semmes, and W. Nazarewicz, Spherical proton emitters, Phys. Rev. C56, 1762-1773 (1997). 



25 



[34] K. Heyde, The Nuclear Shell Model, (Springer Verlag, Berlin, 1990). 

[35] B. Buck, A. C. Merchant and S. M. Perez, Half-lives of favored alpha decays from nuclear ground states, Atomic Data and 

Nuclear Data Tables 54 (1), 53-74 (May, 1993). 
[36] A. D'Arrigo, N. V. Eremin, G. Fazio, G. Giardina, M. G. Glotova, T. V. Klochko, M. Sacchi and A. Taccone, Investigation 

of bremsstrahlung emission in a-decay of heavy nuclei. Physics Letters B332 (1-2), 25-30 (July, 1994). 
[37] V. Yu. Denisov and H. Ikezoe, Alpha-nucleus potential for alpha-decay and sub-barrier fusion, Phys. Rev. C72, 064613, 

p. 9 (2005), nucl-th/0510082 

[38] S. Peltonen, D. S. Delion, and J. Suhonen, a-decay spectroscopy of deformed nuclei reexamined, Phys. Rev. C78, 034608 
(2008). 

[39] T. Ohtsuki, H. Yuki, K. Hirose, T. Mitsugashira, Status of the electron accelerator for radioanalytical studies at Tohoku 

University, Czech. Journ. Phys. 56, D391-D398 (2006). 
[40] V. Yu. Denisov and A. A. Khudenko, Alpha-decay half-lives, alpha-capture and alpha-nucleus potential, Phys. Rev. C79, 

054614 (2009) [23 pages], arXiv:0902.0677 
[41] V. Yu. Denisov and A. A. Khudenko, Alpha decays to ground and excited states of heavy deformed nuclei, Phys. Rev. C80, 

034603 (2009) [10 pages). 



